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O 

, A recently introduced hierarchy of states of a single mode quantised radiation field is examined 

for the case of centered Guassian Wigner distributions. It is found that the onset of squeezing among 
such states signals the transition to the strongly nonclassical regime. Interesting consequences for 
the photon number distribution, and explicit representations for them, are presented. 

> 

I. INTRODUCTION 

Squeezed states of light, and other states exhibiting cither antibunching or subpoissonian photon statistics or both, 
are well known examples of so called "nonclassical" states of radiation jl| Q || . In fact these are the simplest 
and most familiar ones out of an infinite hierarchy of independent signatures of nonclassical states in quantum optics; 
many other signatures have been presented in the literature H . 

The precise definition of a nonclassical state of radiation is based upon the diagonal coherent state expansion of 
Oh, the density matrix p of the state in the quantum theory. Limiting ourselves to the single mode radiation field this 
expansion is |§] 

Jh 

p= I ^(z)\z)(z\, (1.1) 



OO 



FT 

where the coherent states \z) are the familiar normalised eigenstates of the photon annihilation operator a with 
>v>( , complex eigenvalue z and <j){z) is a real normalised weight function which is in general a distribution. The state p is 
5_j ■ said to be "classical " if <j){z) is pointwise nonnegative, and nowhere more singular than a delta function, so that it 
can be interpreted as a classical probability density over the complex plane. Otherwise p is a "nonclassical " state. 
This classification is clearly invariant under rotations and translations in phase space. 

It has been shown elsewhere that there is a dual operator based approach to this distinction between classical and 
nonclassical states, which is physically quite instructive [Q. The representation (1.1), as is well known, is closely 



related to the normal ordering rule of correspondence between classical dynamical variables and quantum operators. 
Given any real classical function f(z*, z) of a complex variable z and its conjugate, one defines a hermitian operator 
F in quantum theory by the replacement z — > a, z* — > and then bringing all factors "by hand" to the left of all 
factors a: 

f(z\ z)^F= f{a\a)\ fit tQ ^ . tQ ^ 

{z\F\z) = f(z\z). (1.2) 
Then the quantum mechanical expectation value of F in the state p is 

(F)=Tr(pF)= f t±<f>( z )f{ z \ z ). (1.3) 
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The key observation now is that while the correspondence / <-> F is linear and takes real functions to hermitian 
operators and vice versa, a real nonnegative f(z*, z) may well lead to a hermitian indefinite F. A state p is then said 
to be classical if this permitted " quantum negativity " in operators never shows up in expectation values, nonclassical 
otherwise: 

p Classical <^ Tr(pF) > for every f(z*, z) > 0, 
p Nonclassical Tr(pF) < for some f(z*, z) > 0. (1.4) 

With this alternative characterization (completely equivalent to the usual one), one has the possibility of defining 
several degrees or levels of nonclassicality, if one restricts in various ways the collection of operators F for which one 
tests the conditions given in (|L4|) @. Specifically, for a single mode system, it has been shown by considering the 
subset of phase invariant (number conserving) operators F which arise from f(z*, z) obeying 

/(z*e^,ze M ) = /(zV) (1.5) 

that an exhaustive and mutually exclusive three-fold classification of states is possible. If f(z*, z) obeys ([l[|), then 
for the expectation value of the corresponding F it suffices to use an angle average of <fr(z): 

oo 

[F, a^a] = Tr(pF) = J dIV(I) f(I 1/2 ,I 1/2 ), 

o 

V(I) = J l^/V*). (1.6) 
o 

One can then obtain the following finer classification of all states: 

p Classical <=> <f>(z) > 0, so V{T) > 0, 
p Weakly nonclassical 4=> V(I) > 0, but <j>(z) 2 
p Strongly nonclassical & V(I) £ 0, so tp(z) ^0 (1.7) 

Thus the previous "nonclassical" has been subdivided now into "weakly nonclassical" and "strongly nonclassical" 
states. Upto and including the weakly nonclassical level, V(I) can be treated as a classical probability density for 
intensity, whether or not 4>(z) can be regarded as a probability distribution over the complex plane; in the third 
strongly nonclassical regime, even V(I) ceases to be a probability density. 

The aim of this paper is to illustrate these ideas in the concrete case of states described by Gaussian- Wigner 
distributions on phase space. It is well known that in a wide variety of physical processes the states of radiation that 
are produced are indeed of this type ||. Their description also lends itself to direct analytical treatment. What we 
shall demonstrate is that within this set of states, the onset of squeezing signals an abrupt change from classical to 
the strongly nonclassical regime; thus the weakly nonclassical states do not show up at all in this family! 

The material of this paper is arranged as follows. In Section II we trace the connection between the descriptions of an 
operator via its Weyl weight and its Wigner representative, and the diagonal weight (j>{z). This gives us a clear picture 
of the extent to which <j)(z) can be a singular distribution, and in turn how singular the quantity V{I) can in principle 
be. Section III examines the class of centered Gaussian Wigner distributions. These are fully parametrised by the 
variance or noise matrix which has to be positive semidefinite and also must obey the uncertainty principle. Among 
these states the only two qualitatively different ones are the nonsqueezed and squeezed ones. In the former case, both 
4>(z) and V(I) can be computed explicitly, and as expected they are finite nonnegative normalized functions. This 
is consistent with their being classified as classical states. In contrast, the squeezed states are shown to be strongly 
nonclassical, and one never sees the weakly nonclassical possibility at all. Section IV gives an example of weakly 
nonclassical states which are naturally outside the Gaussian Wigner family, and offers some concluding remarks. 



II. NATURE OF THE DISTRIBUTIONS 4>{z) AND V{I) 

It is useful to begin by recalling the general properties of the diagonal weight 4>{z) and its angular average P(I), 
and by giving an indication of the kinds of singular distributions we must be prepared to encounter [0] . This is best 
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done by viewing the set of all possible density matrices p as a subset of the family of Hilbert-Schmidt (H-S) operators. 
An operator A on Hilbert space is of H-S type if 

Tr(A f A) < oo, (2.1) 

and among H-S operators we have a natural inner product : 

(A, B) = Tr(A^B). (2.2) 

We deal throughout with systems of one degree of freedom, and with the annihilation and creation operators a, a) 
related to hermitian q and p in the standard way: 

a=-^=(q + ip), a) = -^={q-ip) (2.3) 

The unitary phase space displacement operators are defined by and have the following properties: 

D(<7, t) — exp(iaq — irp), — oo < <7, r < oo; 
£>(<7,r)t = £>(<t,t)- 1 =D(-<7,-t); 
Tr (D(a',T')^D(a,T)) = 2-k5{g' - ct)5(t' - r). (2.4) 

Any H-S operator A can be expanded in the form of an operator Fourier integral representation using its "Weyl 
weight" A(a, r) as expansion coefficient Eg ]: 

f f dadr ~ , . 



A(,7,r) = -L (£>(„, r), A), 

V 27T 

Tr(AU) = (A, A) = [ [ da dr\A(a,T)\ 2 . (2.5) 



Thus the H-S property (2T) of ^4 is translated exactly into the L 2 property of A(a, r) over T?. 2 . 

From -A(<t, r) we pass to the Wigner representative or Wigner distribution W(q,p) of the operator A by a double 
Fourier transform at the c-number level |p"l| : 

W{q 1 p)= / / -^-i(cr,r)exp(icrg-irp) (2.6) 
J J (2tt)5 

Here q and p are canonical coordinates over a classical phase space, and in case A is hermitian its Wigner representative 
W(q,p) is real. Now the H-S property for A amounts to W(q,p) being an L 2 function over 1Z 2 : 

Tr(A^A) = (A, A) = 2tt / J dqdp\W(q,p)\ 2 (2.7) 

For density matrices we are also interested in the ordinary trace: 

Tr (A) = V2^4(0,0) = / / dqdpW{q,p) (2.8) 



It is in the passage from A(a,r) or W(q,p) to </>(z) that the distribution character of the latter shows up. From 
the diagonal representation 

A = Jd^ mlz){zl (29) 
where z = -^{% + iy), when we connect up with the previous relations (2.5, 2.6) we get the result: 

J J \/2tt 



2tt 



J J dqdpWiq.py^" 7 ^. (2.10) 
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Thus the most singular kind of <j>(z) is one whose Fourier transform is the increasing Gaussian factor exp \{ct 2 + t 2 ) 
times a square integrable function A(a,r) - this is the worst behaviour that can in principle occur. Conversely for a 
classical state A(a,r) must more than overwhelm this exponential factor and moreover yield a nonnegative 4>(z). 

Let us next see what this situation for <p(z) entails for its angular average V{I). We work directly with the Wigner 
distribution W(q,p) and find after performing the angular integration: 







dadr 



2tt 



e3 (ff +T )jo(V2I(a 2 +r 2 )) / / dqdpW '(q,p)e l{Tp - aq) . (2.11) 



If we substitute a = V 2K cos tp , r = \/2Ksmip, we can carry out one more angular integration and bring V(I) to 
the following form: 

oo 

V{I) = J dKe% .Jo(2VlK) J J dqd P W(q,p) J W2K(q 2 +p 2 ) 
o 

oo oo 2tt 

= J dKe% Jo{2\flK) J dLJ (2VKL) J d x W(V2Lcosx,V2Lsm X ) (2.12) 



Now just as the relation (2.10) between cj>(z) and A(a, r) involved the classical two dimensional Fourier transformation, 



here one is concerned with the single variable Fourier-Bessel transformation over the half- line (0, oo) which states [ [12[ : 

oo 

J dl\f(l)\ 2 <^^ 



oo 

/(-O = / dKg{K)J {2sflK), 
o 

oo 

g{K) = J dIf(I)J (2VlK), 
o 

oo 

dl\f(l)\ 2 = J dK\g(K)\ 2 , 
o 

J dKJ Q {2\fLK)J Q {2y/lK) = 5(1 - L). (2.13) 
o 

This means that the most singular possible behaviour for which can in principle occur is that its Fourier-Bessel 
transform can be the factor e% times a square integrable function of K over the domain (0, oo) , namely the Fourier- 
Bessel transform of the angular average of W(q,p). The factor e= is just the earlier factor e*' r+r - 1 present in 



eq (2.10); and the situation for V(I) is marginally better than for <j>(z) since now only the angular average of <j>{z) is 
involved. 

The use of phase space language in describing operators in quantum mechanics leads naturally to an examination 
of the behaviours of <j)(z) and V(I) under phase space rotations and translations. As is easy to see, their behaviour 
under rotations is simple: 

W'(q,p) = W(qcosa ~ psma,pcoaa + qsma) (j)' (z) = (f>(ze la ) => V' (I) = V(T). (2.14) 
This invariance of V(I) is as expected. Under translations we have 

W'(q,p)=W(q-q ,p-po)&cf>'(z)=(l>(z- z ), z = -^=(q + ip ) (2.15) 
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However now V'(T) is not expressible in terms of V{I) alone as phase sensitivity is introduced by a translation. There- 



fore while our threefold classification scheme (1.7) is obviously invariant under phase space rotations, the behaviour 
with respect to translations is much more subtle. 

It is evident that the classical states with both <fi(z) and V(I) nonnegative remain classical under translations. 
However a weakly nonclassical state becomes strongly nonclassical for a suitably chosen translation, as the following 
physical argument shows. At the origin V(0) reduces to (f)(0) as no angular average remains. If a weakly nonclassical 
state is given, its (j)(z) must become effectively negative somewhere in the complex plane. By translating the origin 
to such a point and then computing "P'(O) we see that the resulting state is strongly nonclassical. Following a similar 
argument we also see that we can recover 4>(z) in its entirety by subjecting the initial state to all possible phase space 
displacements zq, <j) ' {z) = <\>{z — zq), and then computing the resulting V (7) and collecting the results. 

We conclude this Section by relating the distribution 'P(I) to the photon number probabilities. Indeed these involve 
a complete independent set of phase insensitive quantities and their expectation values: 

* (z*z) n 

f(z*,z) = e-* * ¥J>-~F=\n)(n\, 
p(n) — Tr(pF) — (n\p\n) 

oo 

/jn 
dlVi^e" 1 — (2.16) 



These p(n)'s always give well defined normalised probabilities for finding various numbers of photons, whether or not 
V(I) is itself a probability density. Formally one can invert the above to get V(I) in terms of p(n), as indeed one 
would expect. If we define the generating function q(K) by 

= E hr Knp{n) (2 - 17) 

71=0 

we see that q(K) converges for all real K and is related to V{T) by 

oo 

q(K) - J2 —r Kn / dlV{l)e- 1 - 

n=0 J 

OO 

= J dir^e- 1 J (2VTk). (2.18) 



Using the formula ( 2.15 ) of the Fourier Bessel transformation again we get the inversion 

oo 

V{I) = e 1 J dKq{K)J (2VlK). (2.19) 



In the classical and weakly nonclassical cases, then, the generating function q(K) is itself well behaved and leads 
to nonnegative V(I), but in the strongly nonclassical case, it causes V{I) to be a distribution, or at any rate not a 
probability. 



III. THE CASE OF GAUSSIAN WIGNER DISTRIBUTIONS 



We consider the family of centered Gaussian Wigner distributions, namely those which have vanishing means for q 
and p |13|. The most general such distribution is determined by a real symmetric 2x2 matrix G 



W G (q,p) 



G 



VdetG 



-exp 



A B 
B C 



(q P)G 



(3.1) 
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The condition that W G (q,p) represent a physically realisable quantum mechanical state imposes the following restric- 
tions on G corresponding respectively to normalisability and the uncertainty principle H] : 



G > 0, ie A + C > 0, A = dctG = AC - B 2 > 0; 



G- 1 



1 

-1 



> , ie A + C>0 , A > A 



Combining these we have the complete set of restrictions on G given by 

A + C>0, < A < 1. 
The noise or variance matrix V is defined and given by 



(3.2a) 
(3.2b) 

(3.3) 



V 

(A<z) 2 



(Aq) 2 A(q,p) 
A(q,p) (Ap) 2 



_ 2 G " 2A 



C -B 
-B A 



J J dqdpq 2 W G (q,p), 
A (q,p) = J J dq dp qpW G (q,p), 
(Ap) 2 = J Jdqdpp 2 W G (q,p). 



(3.4) 



Here the vanishing of the means of q and p has been used. In terms of V, the uncertainty principle appears in the 
following form fla]: 



1 1 

detV= — > -. 
4A ~ 4 



(3.5) 



We can use the covariance of <j>(z) and the invariance of V(I) under phase space rotations to simplify the situation 
and to assume without loss of ge ner ality that G and V are diagonal. Moreover these rotations do not disturb the 
three- fold classification of states (1.7). Therefore we parametrise G and V using two real positive parameters a and 
j3 as follows: 



V = 
G = 



1 / a 2 



(3 2 
I /a 2 

i/p 2 



1 



exp 



q_ 

a 2 



, a, j3 > 0, a/9 > 1; 

p 2 



(3.6) 



To deal with <j>{z) and 7 7 (/) we need respectively the Fourier transform and the angular average of W( a m(q,p); these 
are: 



// 



dqdpW( a m(q, p)exp(irp — icrq) = exp 



6 2 



ZTT 

/ dxW^miV^L cosx, VZLsinx) = — jexp ( + -j) ] I (i(— - 

J ap \ a p J \ a p 



(3.7a) 
(3.7b) 



Here Iq(w) = Jo(iw) is the Bessel function of order zero and imaginary argument. 

Returning to the Wigner function W^ a ^{q,p), the nonsqueezed case corresponds to both a, (3 > 1; while if one 
of them becomes less than unity we have a squeezed state. For definiteness in the latter case we take p to be the 
squeezed variable, so we take f3 < 1 and a > 1 maintaining a{3 > 1. Formally we have throughout, on combining 
eqs (pT0|, |3~7aT): 



<A(a,/3)(z) 



dadr 
2tt 



-i( a 2 -l)a 2 -i(/3 2 -l)r 2 



(3.8) 



G 



In the nonsqueezed regime these integrals can be computed and we get expected results: 



2(a 2 - l)" 1 / 2 .^ 2 - 1)-V2 



exp 



y 

p 2 - 1 



a,(3 > 1; 



</»(a,/3)(«) = < 



27T(5(a;) % /2(/3 2 - l)" 1 ^^ 



)-i/2 c 
a = /? = 1. 



27T(5(y)V2(Q 2 - l)- 1 /-exp 



P^l 



a z - 1 



a = > 1; 
a > l,/3 = 1; 



(3.9) 



In all these cases the state is classical. However, once P dips below unity, we see from eq ( |3.8| ) that the Fourier 
transform of <f)(z) is an increasing Gaussian in the variable r. This means that <j)(z) has switched abruptly to being a 
distribution, essentially of the most singular kind that can arise. (Of course, if /? continually decreases and squeezing 
increases, <fi(z) does become more and more singular). This is consistent with squeezed states being nonclassical. The 
interesting point is that there is no intermediate regime (among " Gaussian- Wigner" states) in which the singularity 
of <j>(z) is somewhat milder, say involving finite order derivatives of delta functions. 

To follow the behaviour of Vr a ^{I) as we pass from the nonsqueezed state to the squeezed regime, and when (3 < 1 
to discriminate between the weakly nonclassical and the strongly nonclassical possibilities, we begin by combining 
eqs ( p. 12 3.7b) to get a formal integral expression for V( a ,p){I)'- 



V (a , m (I) = ^ J dKe K ' 2 .M2VlK) J dL e J (2VlK)I q - 4)) 



(3.10) 



The first integral, over L, always converges thanks to the asymptotic behaviours of Jq(z) and Iq(z): 

~2 



-oo \j2/kz 



(3.11) 



Moreover, by suitable and permitted analytic continuation of a standard definite integral available in the litera- 
ture [reference 16 p. 711, formula 6.644] we obtain a formula with whose help the L-integral can be done explicitly. 
The requisite formula is, for real parameters a, b, c obeying a > |c| > 0, b > 0: 



Taking a 



dx e~ ax J {2Vbx) I (cx) 



exp 



-ab 



2 2 

a — c 



Io 



cb 



2 2 

a — c 



(3.12) 



a' 13 



}p,b = K,c= -\ — i here and using the result in eq ( 3.10 ) we get for P( a ,fi) (-0 the single integral 



V [a , p) {I) = JdK e K l\ J (2Vm) e- K{ ^^ ] I (|(a 2 - /3 2 )) 



(3.13) 



First let us look at the classical nonsqueezed situation. Leaving aside the marginal cases when a or f3 equals unity, 
we again use the result (3.12) to evaluate ( 3.1 3| ) explicitly: 



a,/3>l : ^ /3) (/) = 2(a 2 -l)- 1 / 2 .(/3 2 -l)- 1 / 2 



exp 



a - I (3-1 



la 



I 



Q — 1 ~ 1 



(3.14) 



This is explicitly nonnegative, and is consistent with the state being classical. In this case, we can go further and 
obtain a closed form expression for the photon-number probabilities pr a m(n). We have: 



r 



7 



1 

n 

a = 1 
b = 



! V(" 2 - 1)(/? 2 - 1) 
1 1 



= / dle- al l n I (bl), 

1 ) J 

a 2 [3 2 - 1 



a 2 - 1 



rf 2 



1 



{a 2 -l)(fi 2 -lY 



(p 2 



(3.15) 



The resulting integral is a known one leading to an expression in terms of the hypergeometric function [reference 16, 
p.711, formula 6.621] 



dx e~ ax x n I (bx) 



,71+1 " 



n 1 n b 
2 + 2'2 +1;1; 7 



(3.16) 



so the probabilities P( a ^){n) are: 



P{ot,p){n) 



72775 



(a 2 -l)(/? 2 -l) 



2 fl 2 



a^/T - 1 



a,/3 > 1. 



n+l 



1 n 

2' 2 



^77 + 77, tt + 1; i;^ 



The combination 2 of a and /3 does not exceed unity as we have a, j3 > 1: 

1 - z = (a 4 - l)(/3 4 - l)/(a 2 /3 2 - l) 2 . 



(3.17) 



(3.18) 



It is interesting to note that the result (3.17) for P( a ,{3){ n ) is a manifestly nonnegative closed-form expression; in this 
respect it may be contrasted with the expression given earlier in the literature || 

Next let us consider the squeezed regime (3 < l,a > 1/(3. Then the exponential factor e K l 2 in the integral in 
eq (3.13) overpowers the remaining factors: 



e K/2 e -K( a > + 0>)/4 Iq ( K(a 2 _ 



1 



K 



a K{l-p 2 )/2 



ttK 



(3.19) 



This means that Vf a ,m (-0 is n0 longer the Fourier-Bessel transform of a square integrable function of K; it has switched 
abruptly from being a classical probability density for intensity to being a distribution, essentially as singular as is 
permitted by the general considerations of the previous Section! 

There is thus no regime in which V( a .p)(I) remains "classical" while <fi is not - the weakly nonclassical possibility 
is not realised at all in the family of Gaussian- Wigner states. Even though / P( Qlj 3)(/) is a distribution in t he sq ueezed 
regime, we can obtain the photon number probabilities by analytic continuation starting from the result ( 3.17 ) in the 
nonsqueezed case. The justification is the following. At the level of Wigner distributions we know that the probability 
P{a,f3){ n ) is the phase space integral of the product of W( a ^)(q 7 p) and the Wigner function W^ n '(q,p) for the nth 
state of the harmonic oscillator fl7|] : 



p = \n)(n\ =* WW(q,p) = U^e~^ + P*) L n (2(q 2 + p 2 )) ; 
p (aj/9) (n) = 27T f f dqdpW ia , 0) (q,p)W {n Hq,p). 



(3.20) 



Here L n (.) is the nth order Laguerre polynomial. Using the rotational invariance of W^(q,p) and eq. (3.7b) for the 
angular average of W( at p)(q,p), we can reduce P( a ,(3)( n ) to a single radial phase space integral: 



(-1)™ 2 

P(a,p) W = ~ -^7j 27r / dLexp \ —2L - L 



1 



L, HZ ) /,, ( — -^) 



(3.21) 
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This is valid for all a and (3 subject to the standard restrictions a, /3 > 1 a[3 > 1. Since we have symmetry in a and 
/3, we may assume with no loss of generality that a > (3. Then the asymptotic behaviour ( [3.11 ) for Iq(z) as z — > oo 
shows that for large L the integrand here behaves like 



i"~ 1/2 exp{-2L(l + l/a 2 )} 



(3.22) 



Thus the integral ( 3.21 ) is absolutely convergent for all a and /3, and is in fact analytic in these variables (in the 
appropriate regions of the complex planes). 

Having established this, we may now go back to the clo sed e xpression ( 3.17 ) valid in the nonsqueczed case and 
analytically continue it to f3 < 1, a/3 > 1. Now from eq. ( 3.1 8j ) we see that the argument z of the hypergeometric 
function exceeds unity, which lies outside the domain of convergence of the power series expansion of F( 1! ^-, %+l;l;z). 
By analytically conti nuing to z > 1, and keeping track of phases generated in switching from (J3 2 — 1) to (1 — /3 2 ) in 
the prefactors in eq. ( 3.17| ), we find that in the squeezed regime we have different expressions for P( a ,f3){ n ) for even n 
and for odd n: 



P{ a ,(3){n) 



2 [( a 2 -l)(l-/? 2 )]" +1 / 2 1 J r(m + , 1/2) F(m + l/2,m+l/2;l/2;l), 



{a [3 



1) 



71+1 



n±l 
Z 2 



a z (3' 



> 1, a > 



^ r(m + , 3/2) ^(m + 3/2,^ + 3/2; 3/2; | 
13 < 1. 



n = 2m, 
n = 2m -) 



(3.23) 



Once again we have manifestly nonnegative closed form expressions B 

The actual expressions for the first few probabilities show the general trend. We find after simplification that, as 
expected both eq. ( |3.17 ) and eq. ( |3.23j ) give identical functions of a and (3: 

2x{( a 2 + l)(/3 2 + l)}- 1/2 
2(a 2 /3 2 -l) x {(a 2 + l)(/? 2 + l)}- 3/2 
{(a 2 - /3 2 ) 2 + 2(a 2 /? 2 - l) 2 } x {(a 2 + l){(3 2 + 1)}~ 5/2 
(a 2 /? 2 - 1) {3(a 2 - (3 2 f + 2(a 2 /? 2 - l) 2 } x {(a 2 + l)(fi 2 + 1)}~ 7/2 



P{a,l3 


(0) 


V{a.fi 


(1) 




(2) 


V(a,0 


(3) 


V{a.fi 


(4) 




(5) 



{3(a 2 - (3 2 f + 24(a 2 - /3 2 ) 2 (a 2 /3 2 - l) 2 + 8(a 2 /3 2 - l) 4 } x {(a 2 + l)(/3 2 + 1)} 



-9/2 



(a 2 /3 2 - 1) {I5(a 2 - (3 2 f + 40(« 2 - 2 f{a 2 f3 2 - l) 2 + 8(a 2 (3 2 - l) 4 } x {(a 2 + l)(f3 2 + 1)} 



-11/2 



(3.24) 



The appearance of the "uncertainty principle factor" (a 2 /3 2 — 1) in P( a ,f3)( n ) f° r °dd n alone is immediately under- 
standable: when the uncertainty limit is saturated and a(3 — 1 , the Gaussian Wigner function VJ^( Qj i/ a ) (q,p) describes 
the squeezed vacuum, for which it is well known that P( a .i/a)( n ) va nishes when n is odd | ft8| . Conversely, even in 
the nonsqueezed regime, despite the uniform looking expression ( 3.17 ), there is a discrimination between the cases of 
even and odd n which is seen when the hypergeometric function is worked out in detail. In the limit a = j3 = 1, we 
have of course just the vacuum s tate, and then P(i i i)(n) vanishes for all n > 1. This case be seen quite explicitly in 
the expressions displayed in eq. (3.24). 



IV. CONCLUDING REMARKS 



We have examined the class of Gaussian- Wigner distributions for a single mode radiation field in quantum optics 
from the point of view of a recently introduced classification of quantum states into three mutually exclusive types - 
classical, weakly nonclassical and strongly nonclassical. We have found that only the first and third possibilities arise 
in this case, corresponding respectively to the nonsqueezed and squeezed situations. As shown elsewhere, there is an 
interesting class of pure states which give physical examples of the weakly nonclassical type. These are superpositions 
of the number states of the following general type: 

m=e-i^±^=e^)\n), (4.1) 
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where a is any complex number and (3(n) is a nonlinear function of n. Here the photon number probabilities are 
independent of f3(n) and follow the Poisson distribution, so V^(I) is a delta function: 



TV (I) = 5(1- a* a) 



(4.2) 



However on the basis of Hudson's Theorem [|9| it turns out that the Wigner function W^(q,p), which is not Gaussian, 
must be negative somewhere, so in turn <j>(z) cannot be nonnegative. This shows that the states ( |4.l| ) are weakly 
nonclassical. 

Our result that the centered Gaussian- Wigner distributions are never weakly nonclassical has an important physical 
consequence. In the regime a > 1, [3 < 1 which corresponds to quadrature squeezing, since V(I) is not nonnegative the 
nonclassical nature of the state must already show up in properties of the photon number distribution probabilities 
P(a,/3)( n )j i e -j y i a phase insensitive quantities. The simplest such signal, namely subpoissonian statistics, does not 
however display the nonclassicality of the state ||. We find after simple algebra that the Mandel Q-parameter is 
always nonnegative: 



Q{a,0) 



(a^a 2 ) - (qtq) 2 
(fit a) 

= 2 { (a 2 - l) 2 + (/3 2 - l) 2 } /(a 2 + /3 2 - 2) 2 > 



(4.3) 



There are however (infinitely many) other signatures of a nonclassical photon number distribution, some of which are 
local in that they involve only a few contiguous probabilities p(n) . For example we have the result M : 



v(i) > o 



l(n) 



(n 
n = 



- l)p(n 
1,2,3,. 



l)p(n + l) -n{p{{n)f > 0, 



(4.4) 



Therefore if any l(n) is negative for some given state, that is evidence for the strongly nonclassical nature of that 
state. For the states Wt a pjfap), taking a = 2, ^ < /? < 1 as an example, we do find explicitly as shown in 
Figure 1 that 1(2), Z(4), /(6) • • • are negative for some range of values of j3 before turning positive as /3 increases; while 
1(3), 1(5) • • • do not display such nonclassical behaviour. 



0.004 



0.002 



-0.002 



-0.004 




0.5 



0.75 



P 



FIG. 1. Violation of the local conditions on photon number distribution in the squeezed regime 
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It is expected that our conclusions will not be altered drastically if we consider general noncentred Gaussian- Wigner 
distributions. This aspect and other examples of states and the cases of two or more modes, will be taken up elsewhere. 
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